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We study squeezed states of phonons, which allow a reduction in the quantum fluctuations of the 
atomic displacements to below the zero-point quantum noise level of coherent phonon states. We 
investigate the generation of squeezed phonon states using a second order Raman scattering process. 
We calculate the expectation values and fluctuations of both the atomic displacement and the 
lattice amplitude operators, as well as the effects of the phonon squeezed states on macroscopically 
measurable quantities, such as changes in the dielectric constant. These results are compared with 
recent experiments. 
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Introduction. — Non-classical photon states such as var- 
ious forms of squeezed states have attracted much atten- 
tion during the past decade [Q. These novel states are 
attractive because they have new statistical and quantum 
mechanical properties. For instance, some of these states 
can achieve lower quantum noise than the zero-point fluc- 
tuations of the vacuum or coherent states. Thus they 
provide a way of manipulating quantum fluctuations and 
have a promising future in different applications ranging 
from optical communications to gravitational wave detec- 
tion In recent years, squeezed states are also being 
explored in a variety of non-quantum-optics systems, in- 
cluding ion- motion and classical squeezing ||] , molecular 
vibrations [0, polaritons 0,^, and phonons in crystals 
[^"H . Ref. propose a second-order Raman scattering 
(SORS) process for phonon squeezing: if the two incident 
light beams are in coherent states, the phonons generated 
by the SORS are in a two-mode squeezed state. Here we 
consider both the continuous wave case studied in |^] 
and the impulsive case studied in [§. The experimental 
realization of squeezed phonons |^| via a SORS process 
has brought attention to the subject of squeezed phonon 
states p^ . 

Regarding detection methods, [^|J|] proposed that if 
the first-order Raman scattering is either very weak or 
prohibited, the second-order stimulated Raman scatter- 
ing process can be used to generate two-mode phonon 
quadrature squeezed states. Moreover, squeezed phonons 
could be detected by measuring the intensity of the re- 
flected probe light [0J|]. This method has been used to 
detect phonon amplitudes, since the reflectivity is closely 
related to the atomic displacements in a crystal. The 
same argument applies for the transmitivity. Measuring a 
transmitted probe light pulse, Ref. [|| observed squeezed 
phonons produced by an impulsive SORS. The intensity 
of the CW SORS signal for many materials might be too 
weak to be detected with current techniques, but might 
be accessible in the future. 

The SORS process originates from the quadratic term 



in the polarizability change SP^p. The photon-phonon 
interaction V that leads to the SORS process is V ~ 

-3 Eq Po^^-'^'Q^jQ^.irE.^E^p . Here, E,^ 
and E2f3 are electric field amplitudes along a and /3 
directions with frequencies lui and 0^2 • The second- 
order polarizability tensor P^'~'^'' satisfies -P^'"*'"' = 



p-qj 



Recall that the complex normal 
mode operator Qqj of the phonons is related to the 



phonon creation b^_qj and annihilation bqj operators by 



?qj ~ ^qj ' 



qj. If the incident photon fields are not at- 
tenuated we can treat the optical fields as classical waves, 
and also consider the different pairs of ±q modes as inde- 
pendent, and treat them separately. Thus, for one partic- 
ular pair of ±q modes, the complete Hamiltonian for the 
two phonon modes involved in the SORS process has the 
form 0: H^^H^- {4"! i?iai?2/3}QqQ-q , 

where iJq = 7iWq{6jj6q -I- 6^q6_q} is the free phonon 
Hamiltonian for the modes q and — q, ujq = {toi — W2)/2, 
and the branch labels j and j' have been dropped. 

Here we consider two different cases. The first is when 
the incident photons are in two monochromatic beams 
0; i.e., with electric fields Ej = £j cos{ujjt + 4>j); j = 
1,2. In the second case the incident photons are in an 
ultrashort pulse whose duration is much shorter than the 
phonon period ||^. 

Squeezed Phonons via Continuous Wave SORS. — Let 
us now first consider the continuous wave (CW) case. 
Because the photons are monochromatic, we can take a 
rotating wave approximation and keep only the on- 
resonance terms in the Hamiltonian. The off-resonance 
terms only contribute to virtual processes at higher 
orders. This approximation is appropriate for times much 
longer than the phonon period. The simplified Hamilto- 
nian has the form 



= i^q - Aq {5q6_q 

Q/3 



+ C.C.} , 



(1) 
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where (j)i2 and Aq refer to the overall phase and am- 
plitude, respectively, of the product of the 2nd-order 
polarizability and the incident electric fields. Recall 
that P^l^"^ is real, therefore the phase (j>i2 has no q- 
dependence. It originates solely from the two photon 
modes. The Schrodinger equation for the ±q-mode 

phonons is ihdt |i/'q(*)) = 'Hq^\t) |'0q(i)) > and its time- 
evolution operator can be solved by a transformation into 
the interaction picture. The result can be expressed as 



l^q(O)) , (2) 



where (q = — iAq t e~'^'^^^/h . Notice that the second fac- 
tor in the time-evolution operator is a two-mode quadra- 
ture squeezing operator . 

In the CW case considered here, the amplitude of the 
squeezing factor Cq grows linearly with time. However, 
this initial linear growth will be eventually curbed by sub- 
sequent phonon-phonon scattering and optical pump de- 
pletion. In other words, the expression for the squeezing 
factor Cq is valid for times much larger than one phonon 
period, but much smaller than phonon lifetimes (because 
this treatment considers non-decaying phonons). Indeed, 
if this growth rate is not fast enough compared to the 
phonon decay rate, the squeezing effect may never reveal 
itself in an experiment. In addition, the phase of the 
squeezing factor is determined by the phase difference of 
the two incoming light waves. If the ±q phonon modes 
are initially in a vacuum state or in a coherent state, 
the SORS will drive them into a two-mode quadrature 
squeezed state |^]. 

The time evolution operator of all the phonon mode 
pairs (instead of just one pair of ±q modes) that are 
involved in this SORS process has the form U{t) = 
riq t^q(^) ■ Therefore, as long as the photon depletion 
is negligible, all the phonon modes that are involved 
in a SORS process are driven into two-mode quadra- 
ture squeezed states. In other words, squeezing can be 
achieved in a continuum of phonon modes by a CW stim- 
ulated SORS process. 

Squeezed Phonons via Impulsive SORS. — Recently, an 
impulsive SORS process has been used to experimentally 
generate phonon squeezing . Here we treat the problem 
expressing the time evolution operator of the system in 
terms of a product of the two-mode quadrature squeezing 
operator and the free rotation operators Since the 
incident photons are now in an ultrashort pulse, the com- 
plete Hamiltonian can be solved exactly in the limit when 
the optical field can be represented by a (5-function. Such 
an approximation is usually considered when the optical 
pulse duration is much shorter than the optical phonon 
period, which is experimentally feasible with femtosecond 
laser pulses. The Hamiltonian for the SORS can now be 
written as H = J2q {-^q ~ ^'q^it)QciQ-q} , where Aq 
carries the information on the amplitudes of the incom- 



ing optical fields and the electronic polarizability. Notice 
that the light-phonon coupling strength Aq in the CW 
case has units of energy, while Aq here has units of ?i. To 
further simplify the problem, we assume that only ±q 
modes are involved in the process. Such a simplification 
is possible when the photon depletion and the phonon 
anharmonic interaction are negligible, so that different 
pairs of phonon modes are independent from each other. 
The Hamiltonian is now 



(3) 



and the Schrodinger equation for these two phonon 
modes is ifidt\ipq{t)) = H^ltpqit))- This equation can 
be solved by separating the free oscillator terms and the 
two-phonon creation and annihilation terms. The result- 
ing time-dependent wavefunction is 



tu 

iV'q(i)) = exp <{ } exp 



: exp {Cq*M-q - Q\bl^}\M^-)) ■ (4) 



Here — -iX'^e^^^'i^'^ /h. Hence the effect of the opti- 
cal pulse is clear: it first applies a two-mode quadrature 
squeezing operator on the initial state, then rotates the 
state by changing its phase . The state will then freely 
evolve after t — 0^. This result is consistent with Ref. j|] 
where the time-evolution operator is expressed in terms 
of real phonon normal mode operators , instead of the 
complex ones used in this paper. Notice that, in contrast 
to the CW SORS, the phase of the squeezing factor C 
for the impulsive case is fixed by the intensity of the light 
pulse. 

Macroscopic Implications. — Now that we have ob- 
tained the phonon states for both the CW and pulsed 
cases, let us consider the macroscopic implications of 
these states. An experimentally observable quantity O 
which is related to the atomic displacements in the crys- 
tal can generally be expressed in terms of Q^: O — 

0(0) -I- Eq(50/5Qq)Qq -|- . . . = Oo + d + O2 + • ■ • 

where the first term Oq = 0(0) is the operator O when 
all Qq's vanish. An example of an experimentally observ- 
able quantity O is the change in the crystal dielectric con- 
stant 6e due to the atomic displacement produced by the 
incident electric fields. To first order in Qq, Se = 6ei — 

Here ^ 



d{5e) 



+ (&-q 

d{Se)/dQq. Indeed, 



h 

;q is the phase of dO/dQq 
a widely used method to track the phases of coherent 
phonons in the time-domain jl^ ] is based on the obser- 
vation of the reflectivity (or transmission) modulation 
SR (ST) of the sample, which is linearly related to Se — 
the change in the dielectric constant due to lattice vibra- 
tions. The above equation for Se indicates that we can 
define a generahzed Q lattice amplitude operator [^,1): 



,(±q)= (6q + 6lq) 



. This gen- 



eralized lattice amplitude Ug(±q) = 2i?e{Qqe'*'5} is the 
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underlying microscopic quantity related to an observed 
reflectivity or transmission modulation when the linear 
term in Qqj ^ei, exists. 

Since different pairs of ±q phonon modes are un- 
correlated to one another, the fluctuation of Oi{— 
Sci) can be expressed as ((AJei)^) = J2q^>o (^/^f^q) 
\d{6e) /dQqf {Aul{±q)) . Here the state is \ip{t)) = 
U{t)\iP{0)) = UqUq{t)\^ci{0)) in either the CW or the 
impulsive case. We can again focus on a single pair of 
±q modes. In the CW case, using Eq. (^, the fluctuation 



Au^Jiq)) 



[cw] 



= 2{e-^''<' cos^(r!q(t)+,^i2/2) 

+e2'^-sin2(r!q(t)+ 012/2)}, (5) 



where = |(^q| = Aq t/h, riq(t) = ojqt + 7r/4, and 
hereafter (. . .) denotes an expectation value on squeezed 
states, unless stated otherwise. Therefore, at certain 
times, the fluctuation {Au'^{±q))^'^^'> can be smaller 
than 2, which is the vacuum fluctuation level. Fur- 
thermore, all the pair of phonon modes that are driven 
by the stimulated SORS process share the same fre- 
quency: ujq = (uji — W2)/2. Therefore, all the fluctu- 
ations (AMg(±q))(^"'' evolve with the same uiq. Notice 
that there is no dependence on ^'q in the final expression 
of (Awg(±q))('^™\ and the squeezing factor phase (/'12/2 
has no q-dependence, all the pairs of modes involved 
through the SORS share the same phase in their fluctu- 
ations. Therefore there can be squeezing in the overall 
fluctuation {{ASei)'^)^™\ Furthermore, the phase of this 
overall fluctuation can be adjusted by tuning the phase 
difference of the two incoming light beams. 

In the impulsive case |^,^, if the ±q-mode phonons 
are driven into a squeezed vacuum state, the fluctuation 
in Ug(±q) is 

(Au2(±q))' ^ 2{e-^< cos2 n'^{t) + eK sin' n'^{t)} , (6) 



where r' 



X' Jh, and n' Jt) 



r!q(t) 



Tq. Again, the squeezing will reveal itself through 
oscillations in ([A((5ei)]^(q))' which is proportional to 
(Aug(±q))'. Note that these oscillations are essen- 
tially the same as the ones obtained in the CW case. 
However, now the squeezing factor is time- independent. 
Also, the t = phase Tr/4 ~ r'^ in Eq. is q- 
dependent. Eq.(^) can be rewritten as (Aug(±q))' — 
2 {cosh2rq -f- sinh2rq sin(2ajqt — r^)} . For small r'^, 
this becomes (Au2(±q))' = 2{1 + 2r'^ + 2r'^sin{2ujqt - 
Tq)}. This expression has essentially the same form as 
the one obtained in [|: (Qlit)) = (Q^(0)){1 + 2^1 + 
2^qSin(2(jjqt+(y5q)}. The small phase term (^q is neglected 
in p] when computing transmission changes. The differ- 
ence in phases, r'^ versus ipq, is negligible in the limit 
of very small squeezing factor, and originates from the 
different interaction Hamiltonians used here and in 
The interaction term in M is proportional to w^(±q) with 



^'q = (notice that their Qq is real and based on stand- 
ing wave quantization |[Tl[|). Therefore, the interaction 
Hamiltonian in is (in our notation) V oc Ug(±q) oc 
SQqQ-q + Qq + Q-q ■ Howevcr, the last two terms in 
this expression do not satisfy momentum conservation, 
we thus did not include them and kept only QqQ^q in 



our interaction term (this form is also used by Ref. |gT| ) . 

If the linear perturbation 6ei due to phonons is neg- 
ligible, such as in then the second order correc- 
tion 02{= Se2) must be considered. When the phonon 
states are modulated by a SORS, so that the ±q modes 
are the only ones which are correlated, then Se2 — 
dQ — QqQ-q • Let us first focus on one pair 
of ±q modes in the CW case. In a vacuum state, 
(0|QqQ-q|0) = 1 ; while in a squeezed vacuum state 

|0)sq, sq(0|gqQ^q|0),q - ( Au^ (±q) ) (™) /2 , with the 

right hand side given in Eq. (||). Therefore, the expecta- 
tion value of QqQ-q in a squeezed vacuum state is peri- 
odically smaller than its vacuum state value. Let us now 
include all the phonon modes that contribute to de2. In 
a vacuum state, (0|(5e2|0) — J2q9'^^/ i^Q^dQ-q) . On 
the other hand, in a squeezed vacuum state. 



{Se2) 



d\Se) 
dQqdQ^ 



-(A«^(±q))(--) 



(7) 



Since the phase <j)i2/2 has no q-dependence, con- 
tributions from the phonon modes sharing the same 
frequency add up constructively. It is thus possi- 
ble that (562) is periodically smaller than its vacuum 
state value. Similarly, in the impulsive case, {662)' = 
2^^ J2q QQ QQ — (AMg(±q))' ; However, the phase factor 
in {662)' has a q~dependence through r^, so that all the 
phonon modes with the same LUq do not contribute to 
{662)' synchronously. In the CW SORS and in the very- 
small-rq limit impulsive SORS the phase of the expecta- 
tion value (QqQ-q) does not depend on q; this is crucial 
to the experimental observation of modulations in the 
dielectric constant, because this q-insensitivity leads to 
constructive summations of all the q pairs involved. Also, 
at a van Hove singularity a large number of modes con- 
tribute to Se2 with the same frequency and phase, thus 
their effect is larger and easier to observe Q . 

Squeezed Phonons via a finite-width SORS. — Of 
course, real light pulses are not (5-functions. Therefore, 
we have also considered a SORS pumped by a light pulse 
with a finite width (smaller than the phonon period T) 
instead of a (5-function. For a fixed peak height /, we find 
Jirt that the optimal pulse width T°p* that maximizes 
the squeezing effect satisfies T°p^ « T/4.4. This calcula- 
tion indicates that the experiments used a pulse width 
which is nearby the optimal value (r/4.4 « 300/4.4 fs 
« 68 fs « Tp). The calculation can be summarized 
as follows. First, in the impulsive Hamiltonian we replace 
the ^-function by a Gaussian with its width Tp as a vari- 
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ational parameter. Since now the Hamiltonian is time- 
dependent in the interaction picture, we cannot directly 
integrate the Schrodinger equation. Instead, we use the 
Magnus method to obtain the time evolution operator 
and keep only the dominant first term. This approxi- 
mation is valid when the pulse duration is shorter than 
the phonon period. We then calculate the width T°p^ of 
the Gaussian that maximizes the squeezing factor. For a 
constant peak intensity, a pulse that is too narrow does 
not contain enough photons; while it can be proven that 
a pulse which is too long (i.e., with a width comparable 
to T), attenuates the squeezing effect. 

Phonon Squeezing Mechanism. — What is the mecha- 
nism of phonon squeezing in the SORS processes? For 
the CW case, the Hamiltonian is the same as an optical 
two- mode parametric process p2| , with the low frequency 
interference of the combined photon modes as the pump, 
the two phonon modes as the signal and idler. The fre- 
quencies of these modes satisfy Uq + cij_q = lui — lu2. 
The impulsive case is slightly different. Although the 
Hamiltonian is similar to a parametric process, the en- 
ergy transfer from the photons to the two phonon modes 
is instantaneous. The resulting phonon state is a two- 
mode quadrature squeezed vacuum state. Indeed, a reg- 
ular parametric process pumps energy into the signal and 
idler modes gradually, while the impulsive SORS does 
it suddenly. The correlation between the two phonon 
modes, and thus the squeezing effect, is also introduced 
instantaneously. Notice that this mechanism is reminis- 
cent of the frequency-jump mechanism proposed in 
In the impulsive SORS, the frequency of the phonon 
modes has an "infinite" i5-peak change at t — 0, while the 
frequency-jump mechanism has finite frequency changes, 
and squeezing there can be intensified by repeated fre- 
quency jumps at appropriate times. However, as it has 
been pointed out in a finite frequency jump up im- 
mediately followed by an equal jump down results in no 
squeezing at all. 

In conclusion, we have studied theoretically the gen- 
eration of phonon squeezing using a stimulated SORS 
process. In particular, wc calculated the time evolution 
operators of the phonons in two different cases: when 
the incident photons are in monochromatic continuous 
waves, and when they are in an ultrashort pulse. The 
amplitude of the squeezing factor initially increases with 
time and then saturates in the CW SORS case, while 
it remains constant in the pulsed SORS case. In addi- 
tion, the < = phase of the squeezing factor in the CW 
SORS, 012, can be continuously adjusted by tuning the 
relative phase of the two incoming monochromatic pho- 
ton beams, while for the pulsed SORS the phase (oc Aq) 
of the squeezing factor is determined by the amplitude 
of the incoming light pulse For both cases we calculated 
the quantum fluctuations of a generalized lattice ampli- 
tude operator and the second order contribution to the 
change in dielectric constant, which is measurable. For 



the finite- width impulsive case, we computed the optimal 
pulse width, in terms of the phonon period, that maxi- 
mizes the squeezing effect. 
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